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摘 要 : 在 未 来 全 球 变化 的 背景 下 ,我 国 干旱 和 半 和 了 
临 重 大 挑战 ,而 对 干旱 和 半 干 旱 区 的 植被 进行 动态 监测 和 影响 评估 是 一 项 十 分 必要 的 工 
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F 旱 区 域 植被 更 加 敏感 ,平衡 经 济 发 展 和 环境 保 


作 。 通 过 分 析 2000 一 2020 年 浑 善 达 克 沙 地 归 一 化 植被 指数 (Normalized difference vegetation index, 
NDVI) 时 空 变化 特征 ,识别 导致 研究 区 NDVI 变 化 的 驱动 机 制 和 主要 的 驱动 因子 ,采用 情景 分 析 探 
讨 其 植被 未 来 变化 的 轨迹 。 结 果 表 明 : 过 去 20a 间 , 浑 善 达 克 沙 地 NDVI 呈 现 波 动 上 升 趋势 ;放牧 
是 影响 此 地 区 植被 变化 的 主要 因子 ,但 放牧 对 其 NDVI 变 化 的 影响 力 有 逐渐 减弱 的 趋势 ;情景 模拟 
显示 , 浑 善 达 克 沙 地 虽然 向 恢复 方向 发 展 ,但 因 其 脆弱 的 生态 环境 ,使 得 其 植被 仍 受 到 人 类 活动 和 


气候 变化 的 胁 迪 。 研 究 结果 为 浑 善 达 克 沙 地 的 生 


根据 当地 的 实际 情况 提出 相应 的 对 策 与 建议 。 


态 恢复 和 生态 建设 工作 提供 一 定 的 理论 支持 ,并 
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文章 编号 : 


干旱 区 履 盖 超过 40% 的 陆地 表面 ,承载 全 球 近 
40% 的 人 口 趾 ,也 是 对 气候 变化 和 人 类 活动 最 敏感 的 
地 区 之 一 外。 在 可 预见 的 未 来 ,极端 干旱 将 更 加 频 
繁 外 , 干 旱 区 面积 持续 扩展 中。 浑 善 达 克 沙 地 为 距离 
京 津 最 近 的 风沙 源 地 , 自 20 世 纪 80 年 代 末 至 2000 
年 ,其 沙漠 化 快速 发 展 “", 表 现 为 植被 覆盖 降低 、 固 
定 沙丘 活化 等 。2000 年 以 后 , 因 京 津 风 沙 源 等 生 
态 修复 工程 的 实施 ,沙漠 化 的 趋势 得 到 缓解 B 
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素 和 人 类 活动 在 植被 变化 的 主导 作用 。 认 为 降水 
tg 气温 是 控制 植被 分 布 格局 的 主要 气候 因子 ”; 但 
也 有 学 者 认为 人 口 "、 载 冀 量 ”耕地 数量 等 人 类 活 
动因 素 对 植被 的 影响 可 能 大 于 气候 因子 "。 而 事实 
上 ,驱动 因子 间 存 在 着 复杂 的 非 线 性 的 交互 作用 ， 
地 理 探测 器 则 可 通过 探测 空间 分 异性 ”, 量 化 和 识 
别 驱 动因 子 之 间 的 复杂 作用 。 如 孟 琪 等 “应 用 该 
模型 发 现 人 类 活动 与 降水 量 ,坡度 等 共同 影响 此 地 


过 牧 ” 开展 ”等 不 合理 的 土地 利用 ,导致 部 分 地 区 
沙漠 化 加 剧 ” 生态 脆 弱 度 增加 …” ,平衡 经 济 发 展 
和 环境 保护 仍 存 在 挑战 ”。 

对 浑 善 达 克 沙 地 植被 变化 的 研究 集中 在 变化 
趋势 及 其 驱动 机 制 方面 。 前 者 多 采用 归 一 化 植被 
指数 (Normalized difference vegetation index, NDVI)/ 
植被 覆盖 度 '“ \ 净 初级 生产 力 “、 叶 面积 指数 “地 
表 反 照 率 "等 分 析 植 被 覆盖 变化 ;对 于 后 者 , 则 采用 
相关 分 析 “、 回 归 分 析 '“ 和 趋势 分 析 等 ,分 析 植 被 
变化 趋势 与 驱动 因素 的 相关 性 ,并 试图 识别 自然 因 
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区 植被 ; 马 永 桃 等 所 认为 相对 湿度 和 人 类 活动 降 
水 量 和 牲畜 密度 以 双 协 同 为 主 , 且 对 植被 变化 解释 
力 超 过 60%。 

浑 善 达 克 沙 地 生态 环境 脆弱 ,更 易 受 到 气候 变 
化 的 影响 ,在 全 球 未 来 升温 1.5 CC 的 情景 下 (概率 
超过 40%)2 ,加 之 大 规模 的 人 类 活动 干扰 ,将 增加 
此 地 区 植被 未 来 变化 的 不 确定 性 。 为 此 ,本 文 试图 
研究 以 下 问题 : (1) 2000—2020 年 浑 善 达 克 沙 地 
NDVI 的 时 空 变化 趋势 ,(2) 量化 气候 变化 和 人 类 活 
动 在 此 地 区 植被 变化 过 程 中 的 贡献 ,辨识 各 类 因子 
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驱动 强度 是 否 发 生变 化 ,(3) 预测 该 区 植被 覆盖 未 
来 可 能 的 变化 轨迹 。 通 过 研究 ,以 期 为 浑 善 达 克 沙 
地 生态 管理 和 生态 建设 提供 决策 支持 ,以 应 对 未 来 
气候 变化 所 带 来 的 不 利 影响 。 


1 材料 和 方法 


1.1 研究 区 概况 

研究 区 位 于 内 蒙古 中 部 锡林郭勒 草原 最 南端 ， 
地 理 位 置 介 于 112.0°~117.5°E,42.0°~43.5°N 之 间 ， 
东西 长 450 km, 平 均 海拔 1200 m; 属 中 温带 大 陆 性 
气候 ,年 均 气 温 1.5 ;年 降水 量 为 150~450 mm, 东 
西 跨度 大 ,年 降水 量 的 80% 集 中 在 夏季 ;主要 生境 
类 型 包括 固定 沙 地 、 半 固定 沙 地 ,流动 沙 地 ;地 带 
性 土壤 为 栗 钙 土 , 非 地 带 性 土壤 为 风沙 土 ;水 分 条 
件 的 差异 导致 东部 为 稀 树 草原 景观 ,西部 为 荡 漠 
草原 景观 ,植被 逐渐 从 中 生 植物 群落 向 旱 生 植物 
群落 过 渡 ”。 
1.20 数据 来 源 与 处 理 

以 NDVI 代表 浑 善 达 元 沙 地 的 植被 生长 状态 ， 
该 数据 来 源 于 美 家 航空 航天 局 网 站 (https:Wlad- 
sweb.modaps.eosdis.nasa.gov/search ) 的 MOD13Q1 数 
据 产 品 ,采用 HANTS 滤波 Chttps://envi.geoscene.cn/ 
appstore/) 消除 异常 值 的 影响 ,生成 2000 一 2020 年 
NDVI 年 最 大 值 合成 数据 集 。 气 候 数 据 来 源 于 中 国 
气象 数据 网 (http://data.cma.cn) ,在 ArcGIS 10.3 中 采 
用 元 里 金 空 间 插 值 转换 为 栅 格 数 据 (2000 一 2020 
年 ) ;高 程 数据 (DEM) 来 源 于 美国 地 质 调查 局 网 站 
的 SRTM DEM (https://glovis.usgs.org) , 在 ArcGIS 
10.3 中 计算 得 到 坡度 、 坡 向 数据 ;土壤 类 型 数据 来 自 
于 《 址 界 和 谐 土壤 数据 库 V1.2) Chttps://www.fao.org/ 
soils- portal/soil- survey/soil-maps- and- databases/harmo- 
nized-world-soil-database-v12/en ) ;植被 类 型 数据 来 自 
于 《内 蒙古 自治 区 植被 类 型 图 》1:1500000 植 被 类 型 
图 ;土地 利用 数据 来 自 中 国 科学 院 资源 环境 科学 与 
数据 中 心 (http://www.resdc.cn/) ,2000 一 2020 年 每 5 a 
一 期 ;2000 一 2020 年 人 口 密度 数据 来 源 于 WorldPop 
发 布 的 世界 人 口 密度 地 图 (https://www.worldpop.org/ 
geodata/listing?id=76) 。 牲 畜 数 量 由 2000—2020 年 
《锡林郭勒 盟 统计 年 鉴 ) 中 获取 ,采用 基于 NDVI 的 
空间 插值 的 方法 空间 化 ,计算 公式 如 下 : 

SUaggr， 

NDVIsum, 


SU, = NDVI, x (1) 
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式 中 :SU; 为 第 ;年 牲畜 饲养 空间 化 图 层 ; NDVI 为 第 i 
年 NDVI 图 层 ;SUaggr; 为 第 i 年 各 ( 弃 ) 县 牲畜 总 数 图 
层 , 由 牲畜 饲养 统计 数据 与 ( 旗 ) 县 界线 图 层 连接 得 
到 ;NDVIsum, 为 第 ;年 NDVI( 旗 ) 县 区 域 累 加 和 值 图 
层 。 上 述 数据 均 通 过 投影 转换 .裁剪 ,生成 分 辨 率 
23250 m ASHE o 
13 研究 方法 

采用 Theil-Sen 中 值 趋势 (Theil-Sen median trend 
analysis , T-S 4} Hf ) fll M-K $% 3m (Mann- Kendall test) 
分 析 浑 善 达 克 沙 地 NDVI 时 空 变化 特征 ,使 用 地 理 
探测 需 识 别 影响 该 沙 地 NDVI 变 化 的 人 为 和 自然 因 
素 ,结合 盘 选 出 的 驱动 因子 ,利用 CA-Markov 模型 对 
该 区 NDVI 进 行 多 情景 预测 。 
1.3.1 Theil-Sen 中 值 趋势 分 析 T-S 分 析 是 一 种 非 
参数 统计 方法 ,适用 于 存在 异常 值 的 时 间 序 列 趋势 
分 析 ” ,公式 如 下 : 


= 
S= medio : 


X; 
a 1<i<j<n (2) 


式 中 :5 为 时 间 序 列 数据 的 斜率 ,反映 数据 的 变化 趋 
势 , 当 S>0 时 ,NDVI 呈 上 升 趋势 变 化 ,S<0 时 ,NDVI 
旺 下 降 趋势 变化 ;median 为 时 间 序 列 的 中 位 数 函 
数 ;x 和 %% 分 别 为 时 间 序 列 i; 年 和 j 年 的 值 ;n 为 此 研 
究 中 时 间 序 列 数据 的 长 度 。M-K 检 验 用 于 检验 序列 
的 变化 趋势 和 突变 出 现时 间 ””"。 选 取 检 验 显著 性 
水 平 a=0.05 ,结果 划分 为 显著 变化 (IZI>1.64) 和 不 显 
著 变 化 (IZI<1.64)。 
1.3.2 驱动 因子 提取 考虑 到 数据 的 可 获取 性 、 区 
域 代表 性 ,并 参考 相关 文献 3, 从 自然 因素 和 人 类 
活动 等 方面 选取 11 种 检测 因子 ( 表 1), 并 以 每 5a 为 
1 个 研究 时 期 ,探测 影响 研究 区 NDVI 分 布 及 变化 的 
影响 因子 。 使 用 ArcGIS 10.3 ff) Fishnet 工具 ,在 区 域 
内 生成 2 kmx2 km 网 格 ,提取 网 格 中 心 点 坐标 和 栅 
格 值 , 代 和 人 地 理 探测 器 中 计算 。 
1.3.8 地 理 探测 人 器 选择 地 理 探测 器 (http://www. 
geodetector.cn) 中 两 种 探测 需 , 以 揭示 研究 区 NDVI 
变化 的 驱动 机 制 。 

(1) 因子 探测 :将 NDVI 作 为 因 变 量 , 各 因子 作 
为 自 变 量 , 则 影响 区 域 NDVI 分 异 决定 值 g 为 : 

Saa 


q=1- += Pp BEL Sy SEL (3) 
no 


式 中 :9 为 度量 空间 分 异性 的 指标 沁 为 分 类 数目 ;mm 
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表 1 浑 善 达 克 沙 地 NDVI 影 响 因 子 
Tab.1 Influence factors of NDVI in Otindag Sandy Land 


xm TENET "TNI 

自然 因素 气候 年 均 温 /?C TEM 
年 降水 量 /mm PRE 

年 均 风 速 /ms WND 

地 形 高 程 /m ELEV 
坡度 /(9) SLP 
坡 向 /(9) ASP 
土壤 土壤 类 型 SOIL 
植被 植被 类 型 VGT 
人 类 活动 土地 利用 类 型 LUC 
AH JEA kim POP 

牲畜 数量 / 头 SHEEP 


和 n SPA Je h ME KAY EAS POR ; on Flo? 分 别 
为 h 和 全 区 的 方差 。g 的 取 值 范围 为 [0, 1] ,9 值 越 
大 ,说 明 影 响 因子 对 NDVI 变 化 的 解释 力 越 强 ?1。 
(2) 交互 作用 探测 :用 于 识别 驱动 因子 之 间 的 
交互 作用 。 通 过 计算 两 种 因子 的 g 值 之 和 与 影响 
子 两 两 联合 的 g 值 ,来 判断 2 个 因子 的 交互 作用 对 
NDVI 的 影响 是 增强 还 是 减弱 ,或 者 2 个 因子 是 独立 
起 作用 (图 1)。 
1.3.4 植被 变化 趋势 模拟 与 多 情景 预测 CA-Markov 
模型 既 保 持 了 CA 模拟 复杂 系统 空间 的 能 力 又 发 挥 
了 Markov 模 型 进行 长 期 预测 的 功能 ”。CA 通常 包 
括 单元 .状态 .邻近 范围 转换 规则 4 个 基本 要 素 , 表 
达 式 如 下 : 


Sai fS N) (4) 

IRP Su SIA H1 ARRERA; 29 

局 部 空间 元 胞 状态 的 转化 规则 ;NN 为 元 胞 的 邻 域 。 

若 一 个 过 程 存在 有 限 多 个 时 序 , 其 任意 时 刻 4, 的 元 

胞 a 只 与 其 前 一 时 刻 的 元 胞 4, 相关, 可 认为 该 过 
程 存在 马尔 科 夫 性 , 则 CA-Markov 模 型 可 以 表示 为 : 

Sa P XS, (5) 


9 交互 作用 : qX1NX2) 


式 中 :Pj 为 状态 转移 矩阵 。 

采用 GeoSOS-FLUS 软件 (http://www.geosimula- 
tion.cn/flus.html) 进 行 多 情景 分 析 。 因 模型 输入 数据 
需 使 用 离散 型 栅 格 , 故 将 NDVI 分 为 5 个 级 别 :高度 
植被 覆盖 区 (NDVI>0.6) 、 较 高 植被 覆盖 区 (0.5S<ND- 
VI<0.6) .中 度 植被 覆盖 区 (0.4<NDVI<0.5) 、 较 低 植 
被 覆盖 区 (0.3<NDVI<0.4) 、 低 度 植 被 覆盖 区 (ND- 
VI<0.3)。 以 2010 年 和 2015 年 的 NDVI 空 间 分 布 图 
作为 基础 数据 ,计算 得 到 NDVI 转 移 和 矩阵 和 适宜 性 
图 谱 以 作为 转换 规则 ,模拟 2020 年 研究 区 的 NDVI 
分 布 ,并 与 真实 值 作 Kappa 精 度 分 析 , 以 检验 模型 在 
植被 变化 动态 预测 中 的 可 靠 性 。 

基于 地 理 探 测 需 的 分 析 结 果 ( 表 3 ,选择 标准 为 
q>0.3) ,利用 CA-Markov 模 型 ,通过 调控 主导 驱动 
子 (自然 因子 和 人 类 活动 因子 ) 范 围 值 的 大 小 ,将 
2030 年 作为 目标 年 , 设 定 以 下 7 种 发 展 情景 ( 表 2)。 


2 结果 与 分 析 


21 浑 善 达 克 沙 地 NDVI 变 化 时 空 分 布 特征 
2000—2020 ^F. NDVI 提高 区 域 的 面积 (77.97%) 
大 于 减少 区 域 的 面积 (11.11%), 其 余 10.92% 的 区 域 
变化 不 明显 。 从 空间 分 布 上 看 , NDVI 提高 区 域 分 
布 在 浑 善 达 克 沙 地 东北 部 、 东 南部 和 中 部 ,NDVI 减 
少 区 域 呈 斑 块 状 分 布 在 中 部 和 东部 地 区 (图 2a)。 
从 时 间 尺 度 看 (图 2b) ,研究 区 NDVI 呈 现 波动 上 升 
趋势 ,最 大 值 出 现在 2018 年 和 2019 年 ,但 其 NDVI 
仍 处 于 较 低 水 平 (均值 为 0.19)。 
2.2 浑 善 达 克 沙 地 NDVI 变 化 驱动 机 制 分 析 
2.2.1 不 同时 期 探测 因子 变化 因子 探测 结果 ( 解 
释 力 g 值 ) 反 映 各 因子 对 浑 善 达 克 沙 地 NDVI 影 响 的 
大 小 。 总 体 来 看 ( 表 3),2001 一 2020 年 各 时 期 g 值 最 
强 的 因子 为 牲畜 数量 (SHEEP) ,平均 值 0.65; 其 他 因 


非 线性 减弱 单 因子 非 线性 减弱 双 因 子 增强 独立 非 线 性 增强 
V 4 V M V 
t t. t 
Min[q(X1), q(X2)] Max[q(X1), q(X2)] q(X1)*q(X2) 
注 :4 为 度量 空间 分 异性 的 指标 ;Z1 1X2 分别 为 导致 NDVI 变 化 的 影响 因子 1 和 2。 


图 1 2 个 变量 对 区 域 NDVI 变 化 的 交互 作用 (修改 自 王 劲 峰 等 ”) 
Fig. 1 Interaction of two variables on NDVI changes (Modified from Wang et al?) 
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表 2 多 情景 预测 的 参数 设置 
Tab.2 Parameters for multi-scenario prediction 
情景 情景 名 称 情景 变化 参数 
惯性 发 展 情景 惯性 发 展 驱动 因子 不 做 处 理 , 遵 循 NDVI 的 惯性 发 展 规律 
气候 变化 情景 ' 降水 增多 调整 降水 值 ,将 降水 的 值 总 体 提升 25% ,其 余 因 子 不 做 处 理 
降水 降低 调整 降水 值 ,将 降水 的 值 总 体 降低 25% ,其 余 因 子 不 做 处 理 
风速 提高 调整 风速 值 ,将 风速 的 值 总 体 提升 25% ,其 余 因 子 不 做 处 理 
风速 降低 调整 风速 值 , 将 风速 的 值 总 体 降 低 25% ,其 余 因 子 不 做 处 理 
经 济 优先 情景 放牧 压力 增加 调整 牲畜 数量 ,将 牲畜 数量 的 值 总 体 提升 50% ,其 余 因子 不 做 处 理 
生态 保护 情景 * 放牧 压力 减少 调整 牲畜 数量 ,将 牲畜 数量 的 值 总 体 降低 50% ,其 余 因 子 不 做 处 理 
注 :1 表 示 根 据 IPCC 第 六 次 报告 汪 和 地 理 探测 器 的 分 析 结 果 ( 表 3) ,并 做 适当 调整 ;2 表示 根据 研究 区 多 年 牲畜 数量 变化 特征 并 做 适当 调整 ;NDVI 为 归 


一 化 植被 指数 。 
(a) NDVI 空间 变化 趋势 (b) 年 均 NDVI 变 化 趋势 
年 均 NDVI 变化 趋势 
026 eene 线性 趋势 
0.24 y=0.0014x+0.1787 
B R?=0.2981 
: E 
图 例 "E ür 
变化 趋势 变化 斜率 () Z 面积 比 /% 
E] 显著 降低 S<-0.0005  Z<-1.64 2.03 一 一 研究 区 边界 
E] 轻微 降低 $«-0.0005 |Z|x1.64 9.08 
变化 不 明显 ， |$|<0.0005 |Z|<1.64 10.92 0 50 km 
国 轻微 提高 S>0.0005  |Z|«1.64 42.68 
B 显著 提高 S>0.0005 Z>1.64 3529 
注 : NDVI 为 归 一 化 植被 指数 ;Z 为 统计 的 值 。 下 同 。 
图 2 2001 一 2020 年 浑 善 达 克 沙 地 NDVI 变 化 趋势 
Fig. 2 NDVI change trend of Otindag Sandy Land from 2001 to 2020 
表 3 2001 一 2020 年 浑 善 达 克 沙 地 探测 因子 gq 值 变化 
Tab.3 Variation of q value of driver factors in Otindag Sandy Land from 2001 to 2020 
HF 
年 份 
TEM PRE WND ELEV SLP ASP SOIL VGT LUC POP SHEEP 
2001—2005 0.09 0.54 0.53 0.49 0.08 0.04 0.25 0.08 0.14 0.05 0.74 
2006—2010 012 0.49 0.50 0.42 0.09 0.04 0.25 0.07 0.14 0.02 0.60 
2011—2015 016 0.54 0.53 0.47 0.08 0.04 0.24 0.08 0.13 0.04 0.64 
2016—2020 020 0.57 0.54 0.49 0.09 0.03 0.25 0.08 0.13 0.08 0.62 


子 排序 为 年 降水 量 (PRE)> 年 均 风 速 (WND)> 高 程 
(ELEV) ,它们 在 各 阶段 o 值 均 在 0.40 以 上 ;土地 利 
用 类 型 (LUC) 年 均 温 (TEM) | EA (SOIL) X} 
NDVI 影 响 一 般 (g 值 平均 为 0.12~0.25); 坡 向 (ASP)、 
人 口 密度 (POP) .植被 类 型 (VGT) 及 坡度 (SLP) 对 
NDVI 影 响 的 4 值 较 低 (平均 为 0.02~0.09)。 

由 各 研究 时 有 段 来 看 ,2001 一 2020 4E TEM , PRE 


和 POP 的 g 值 总 体 呈 上 升 趋势 ,而 SHEEP 的 g 值 呈 
现 降低 的 趋势 (" 值 由 0.74 降 低 到 0.62 ) ,其 他 因子 变 
化 幅度 不 明显 。 意 味 着 近年 来 气候 因子 对 NDVI 的 
影响 力 增 大 ,而 放牧 等 人 类 活动 对 NDVI 的 影响 力 
2.2.2 不 同时 期 探测 因子 交互 作用 不 同 研究 时 段 

下 各 探测 因子 交互 作用 下 的 g 值 均 大 于 任意 单 因 子 
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的 4 值 ( 图 3) ,表现 为 双 因 子 增强 和 非 线 性 增强 ,不 
存在 独立 或 减弱 关系 。 在 2001 一 2005 年 和 2006 一 
2010 年 这 2 个 时 段 ,WND FI SHEEP 的 交互 作用 最 
强 (g 值 分 别 是 0.82 和 0.69); 而 在 2011 一 2015 年 和 
2016 一 2020 年 这 2 个 阶段 ,SOIL 和 SHEEP 的 交互 作 
用 最 强 (g 值 均 为 0.73)。 

在 整个 研究 时 段 中 ,各 因子 间 的 交互 作用 亦 发 
生变 化 ,表现 为 气候 因子 与 其 他 因子 的 作用 增强 、 
人 类 活动 因子 与 其 他 因子 交互 作用 减弱 。 不 同时 
期 下 气候 因子 间 交 互 作用 下 的 g 值 增幅 显著 ,表明 
气候 因子 间 交 互 协同 作用 对 NDVI 影 响 较 大 ,说明 
气候 因子 间 是 相互 作用 而 非 独 立 影响 NDVI 的 空间 
分 布 , 具 体 表 现 为 : 除 SHEFEP 外 ,TEMPRE 与 其 他 
因子 交互 作用 的 g 值 均 有 小 幅 的 增加 (幅度 0.01~ 


(a) 2001 一 2005 年 


TEM PRE WNDELEVSLP ASP SOILVGT LUC POP 
因子 
(c) 2011 一 2015 年 


0.10) ; WND 则 表现 出 与 SOIL、LUC 和 SHEEP 等 因 
子 交 互 作用 减弱 的 趋势 ,但 变化 幅度 较 小 ,分 别 减 
少 0.03 .0.04 和 0.09。 在 人 类 活动 因子 中 ,POP 表现 
出 与 其 他 因子 的 交互 作用 增强 的 变化 趋势 ,但 g 值 
变化 幅度 较 小 ;而 SHEEP 与 其 他 因子 间 的 交互 作用 
4 值 在 整个 研究 时 段 的 有 减少 趋势 (减少 幅度 分 别 
H 0.08-0.13) ,表明 随 着 时 间 的 推移 ,放牧 对 NDVI 
的 影响 逐渐 减弱 ,但 放牧 依然 是 影响 NDVI 变 化 的 
主要 原因 。 地 形 .土壤 和 植被 等 因子 间 交 互 作用 g 
值 变 化 不 明显 (在 +0.01 间 波动 )。 
2.3 基于 CA-Markov 模 型 的 NDVI 时 空 变化 模拟 
及 预测 

通过 对 比 2020 年 实际 NDVI 数 据 与 模拟 NDVI 
数据 (图 4) ,模拟 总 体 精度 为 65% 、 最 佳 优 值 系 数 


(b) 2006—20104£ 


TEM PRE WNDELEV SLP ASP SOIL VGT LUC POP 
因子 
(d) 2016 一 2020 年 


TEM PRE WND ELEVSLP ASP SOIL VGT LUC POP 
因子 


TEM PREWNDELEVSLP ASP SOIL VGT LUC POP 
因子 


4 值 
0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00 


注 :TEM 为 年 均 温 ;PRE 为 年 降水 量 ;WND 为 年 均 风 速 ;ELEV 为 高 程 ;SLP 为 坡度 ;ASP 为 坡 向 ;SOIL 为 土壤 类 型 ;VGT 为 植被 类 型 ; 
LUC 为 土地 利用 类 型 ;POP 为 人 口 密度 ;SHEEP 为 牲畜 数量 。* 表 示 2 因子 为 非 线性 增强 关系 ,其 他 数值 均 表示 2 因子 为 


E 


因子 增强 关系 ,结果 均 通 过 a=0.05 水 平 的 检验 。 


图 3 2001 一 2020 年 浑 善 达 克 沙 地 各 阶段 交互 探测 器 结 
Fig.3 Results of interactive detectors in each stage of Otindag Sandy Land from 2001 to 2020 
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(Figure of merit, FOM) 为 0.13, 模 拟 结 果 较 好 。 各 情 
景 下 NDVI 的 空间 分 布 如 图 5 所 示 : 各 等 级 NDVI 分 
布 趋势 总 体 上 保持 不 变 , 低 度 植被 覆盖 区 集中 分 布 
在 研究 区 西部 沙丘 地 区 ,高 度 植被 覆盖 区 则 位 于 东 
部 地 势 平坦 地 区 ,各 情景 中 不 同 植被 窗 盖 区 相互 转 
化 的 区 域 主要 位 于 研究 区 的 中 西部 。 在 惯性 发 展 


(a) 实际 NDVI 


植被 覆盖 
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低 度 植被 覆盖 区 面积 占 比 最 大 (10527 km’, 

积 占 比 为 23% ) ,其 次 是 中 等 植被 覆盖 区 和 较 低 
ee REE etd ie hi ome 但 较 
2020 年 植被 覆盖 整体 提高 ,表现 为 中 等 植被 覆盖 区 
和 较 高 植被 覆盖 区 面积 呈 增 加 趋势 ,分别 增加 了 
1445 km 和 1309 km:( 表 4)。 在 气候 变化 情景 下 , 降 


(b) 模拟 NDVI 


| 
图 例 0 50km 


园 低 度 (<0.3) 国 较 低 (0.3~0.4) 加 中 等 (0.4~0.5) 国 较 高 (0.5~0.6) B] 高 度 C0.6) 一 一 研究 区 边界 


图 4 2020 年 浑 善 达 克 沙 地 实际 NDVI 与 模拟 NDVI 空 间 分 布 
Fig.4 Spatial distributions of actual NDVI and simulated NDVI in Otindag Sandy Land in 2020 


(a) 惯性 发 展 


nu 


(c) 降水 降低 


(e) 风速 降低 


i 


(g) 生态 保护 


m. di 


(b) 降水 增多 


maj 


(d) 风速 提高 


i 


0 100 km 人 
a = 


B co.) B 较 高 0.5~0.6) 
[eo c.o 
E 中 等 (0.4~0.5) — — 研究 区 边界 


图 5 2030 年 浑 善 达 克 沙 地 各 情景 NDVI 模拟 空间 分 布 
Fig. 5 Spatial distributions of simulated NDVI under different scenarios in Otindag Sandy Land in 2030 
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水 增多 情景 和 风速 降低 情景 的 植被 覆盖 整体 提高 ， 
较 高 植被 覆盖 区 和 高 度 植被 覆盖 区 面积 分 别 增加 
了 1187 km' fil 2855 km 1309 km 和 1311 km’; 而 降 
水 降低 情景 和 风速 提高 情景 下 的 植被 覆盖 整体 降 
AI , 低 度 植被 覆盖 区 面积 分 别 增加 了 271 km? 和 


3)。PRE 是 影响 NDVI 变 化 的 主要 气候 因子 , 它 决定 
着 干旱 区 - 半 干 旱 区 植被 的 空间 分 布 和 生长 状态 ,为 
主要 限制 因子 ,这 与 马 永 桃 等 *、 李 永利 等 ”Sun 
等 "Ma 等 "在 浑 善 达 克 沙 地 的 相关 研究 结论 相 一 
致 。 人 研究 结果 显示 WND 对 研究 区 NDVI 变化 影响 


724 km 。 在 经 济 优 先 情景 下 ,增加 牲畜 饲养 数量 ， 
模拟 较 高 放牧 压力 情景 下 植被 的 演变 状况 ,与 2020 
年 相 比 , 低 度 植 被 覆盖 区 和 较 低 植被 覆盖 区 面积 4 
别 增 加 了 716 kw 和 162 km , 即 放牧 会 导致 植被 覆 
盖 降 低 。 在 生态 保护 情景 下 ,放牧 密度 降低 ,草地 
受到 人 类 活动 的 破坏 减少 , 较 高 植被 覆盖 区 和 高 度 
植被 覆盖 区 面积 分 别 增加 了 1309 km 和 2995 km’, 
植被 覆盖 增加 面积 明显 高 于 降水 增多 情景 和 风速 
降低 情景 ,进一步 说 明 人 类 活动 对 NDVI 的 影响 较 
单一 气候 影响 大 。 


3 讨论 


3.1 浑 羡 达 克 沙 地 NDVI 变 化 原因 

2000 一 2020 年 浑 善 达 克 沙 地 的 NDVI 呈 现 波动 
上 升 的 趋势 ,增加 趋势 为 平均 0.19 a, A NDVI =Œ FA 
西向 东 递 增 趋势 ,这 与 马 永 桃 等 六 、 安妮 等 5 的 研 
究 结果 相似 。20 a 间 , 浑 善 达 克 沙 地 向 恢复 方向 发 
展 ,但 其 脆弱 的 生态 环境 仍 受 到 气候 变化 和 人 类 活 
APSE AY BO 

分 析 发 现 , 浑 善 达 克 沙 地 NDVI 变 化 受到 人 类 
活动 和 自然 因子 的 影响 ,其 中 以 人 类 活动 为 主 ,gq 值 
最 强 的 因子 为 SHEEP( 表 3)。 较 多 研究 均 发 现在 
我 国 北方 草原 区 牲畜 数量 与 NDVI 有 很 强 的 相关 关 


较 大 ,这 可 能 与 其 地 貌 特征 有 关 , 如 陈 宸 等 中 也 发 
现 平 均 风 速 是 新 疆 荒 漠 草 地 NDVI 的 重要 影响 
子 。 在 植被 较 好 的 地 区 (固定 沙 地 和 低 湿 地 等 ) , 风 
速 增加 导致 土壤 莱 散 量 增加 ,增加 植物 的 水 分 胁 
3B ,影响 草地 NDVI” ;在 植被 较 差 的 区 域 ( 半 固定 
和 流动 沙 地 等 ) ,过 高 的 风速 吹 走 土 壤 、 损 伤 植物 幼 
苗 等 信守。 此 外 ,因子 探测 分 析 中 发 现 ,TEM 对 ND- 
VI 解释 能 力 较 弱 (4 个 研究 时 段 g 平 均值 在 0.09~0.20 
之 间 )。 这 与 李 永 利 等 "认为 此 沙 地 地 上 净 初 级 生 
产 力 与 年 极端 最 高 气温 呈 极 显著 相关 关系 (r=-0.4) 
的 结果 相左 。 可 能 的 原因 是 本 文采 用 平均 温度 指 
标 ( 历 年 变化 缓慢 ) , 而 不 像 极 端 气温 增加 /变化 明 
显 , 因 此 导致 TEM 对 NDVI 变 化 的 解释 能 力 较 弱 。 
另外 ,也 有 学 者 认为 除 春季 外 ,中 国 北方 草原 区 生长 
季 各 季节 植被 对 温度 的 响应 不 如 降水 敏感 ““; 王 旭 
洋 等 中 认为 中 国 北方 四 大 沙 地 NDVI 变 化 对 气温 的 
变化 不 敏感 ,但 也 受到 气温 与 降水 的 协同 驱动 中。 
本 研究 也 支持 该 观点 ,TEM 单 因子 对 NDVI 变 化 解 
释 能 力 较 弱 (多 期 9 值 均 值 为 0.14) ,而 TEM 与 PRE 
的 交互 作用 对 NDVI 变 化 影响 较 大 ,各 人 研究 时 段 o 值 
均 在 0.54 以 上 (图 3)。 

较 长 时 间 尺 度 上 ,气候 因素 控制 植被 覆盖 分 布 
格局 ,而 在 短 时 间 内 人 为 因素 对 植被 覆盖 的 影响 比 


系 ,是 导致 植被 变化 的 主要 驱动 因素 在 自 
然 驱 动力 方面 , NDVI 受到 PRE 和 WND 的 影响 ( 表 


气候 因素 更 为 严重 "。 本 研究 发 现 2000 年 后 ,虽然 
放牧 是 影响 研究 区 植被 变化 的 主要 因子 ,但 气候 因 


表 4 浑 善 达 克 沙 地 现状 及 模拟 各 等 级 NDVI 面 积 


Tab. 4 Area of actual and simulated NDVI of Otindag Sandy Land under different scenarios Ikm? 

年 份 情景 情景 名 称 - 植被 覆盖 aia = H 
低 度 较 低 中 等 较 高 高 度 
2020 - - 12034 10105 8531 5607 5469 
2030 ARERR S 惯性 发 展 10527 9141 9976 6916 5187 
(模拟 ) 气候 变化 情景 降水 增多 7511 9141 9976 6794 8324 
降水 降低 12305 10729 7992 5532 5188 
风速 提高 12758 10227 7990 5760 5011 
风速 降低 8933 9141 9976 6916 6780 
经 济 优先 情景 放牧 压力 增加 12750 10267 7927 5499 5303 
生态 保护 情景 放牧 压力 减少 7249 9141 9976 6916 8464 
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子 对 NDVI 的 影响 力 逐 渐 增 大 ,放牧 等 人 类 活动 对 
NDVI 的 影响 力 在 逐渐 减弱 ( 表 3、 图 3)。 这 可 能 与 生 
态 恢复 工程 及 相关 的 实施 草 裔 平 衡 政 策 有 关 ? VU 
2000 年 以 后 ,相继 开展 的 “ 京 津 风 沙 源头 治理 工程 ” 
等 生态 工程 的 ,对 区 域 植被 覆盖 变化 产生 了 重要 的 
影响 站。 且 当 地 政府 实施 全 区 禁 牧 、 休 牧 、 草 畜 平 
衡 工作 ,也 促使 植被 恢复 ”。 
3.2 政策 启示 

浑 善 达 克 沙 地 NDVI 呈 现 “ 东 高 西 低 ” 的 分 布 特 
征 ,植被 改善 的 区 域 主要 集中 在 浑 善 达 克 东 南部 、 
西部 地 区 ,植被 改善 不 明显 。 结 合 情 景 分 析 , 本文 
认为 ,采取 因地制宜 地 分 区 治理 方式 ,对 不 同 的 分 
区 采取 不 同 的 治理 对 策 。 在 植被 较 好 的 东部 地 区 ， 
利用 与 保护 并 举 ,遵循 “以 水 定 草 We” B JER 
则 ,科学 布局 草 业 和 畜牧 业 。 把 解决 家 畜 超 载 问题 
作为 遏制 沙 地 退化 的 首要 措施 ,根据 降水 的 变动 情 
况 ,动态 判断 草场 承载 能 力 。 同 时 ,把 发 展 畜 牧 业 
产业 化 作为 农 牧 业 产业 化 突破 口 , 有 条 件 地 区 实施 
人 工种 草 推进 畜牧 产业 化 。 或 通过 产业 结构 调整 
LAE SIRE UE BI Az REO ,发挥 区 位 优势 ,积极 发 
展 沙 地 旅游 . 康 养 产业 等 。 对 于 植被 欠 佳 的 中 西部 
地 区 ,以 保护 为 主 利用 为 辅 。 可 借鉴 牧民 的 传统 放 
牧 方式 ,将 沙 地 作为 冬季 营地 使 用 可 有 效 保护 沙 地 
植被 ”; 或 实行 跨 区 域 草 畜 平衡 ,配合 周边 的 典型 
草原 区 实行 季节 性 倒 场 放牧 或 从 周边 农 区 调 草 ; 充 
分 发 挥 “ 京 津 综 ”生态 服务 供应 地 的 作用 ,探索 基于 
生态 系统 服务 的 转移 文 付 ,减少 放牧 对 植被 的 影 
响 , 以 应 对 未 来 气候 变化 所 造成 的 不 利 影响 。 


4 结论 

(1) 2000—2020 年 浑 善 达 克 沙 地 的 NDVI 呈 波 
动 上 升 的 趋势 ,年 增加 0.19。NDVI 变 化 呈 空 间 异 质 
性,35.29% 的 区 域 呈现 显著 改善 的 趋势 ,显著 改善 
的 大 部 分 地 区 集中 在 浑 善 达 克 沙 地 东南 部 ,地 表 植 
被 改善 的 面积 大 于 植被 退化 的 面积 ,植被 显著 退化 
的 地 区 主要 分 布 在 研究 区 的 东部 。 

(2) 放牧 仍 是 影响 研究 区 NDVI 变化 的 主要 因 
子 , 但 气候 因子 对 NDVI 的 影响 力 增 大 ,放牧 对 ND- 
VI 的 影响 力 在 逐渐 减弱 。 

(3) 虽然 浑 善 达 克 沙 地 向 恢复 方向 发 展 ,但 其 
脆弱 的 生态 环境 仍 亦 受到 人 类 活动 和 气候 变化 的 
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胁迫 。 多 情景 模拟 结果 为 制定 统筹 区 域 社会 经 济 
发 展 与 生态 治理 的 草地 恢复 政策 提供 借鉴 和 数据 
xt 
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Vegetation dynamic and its driving force and multi-scenario 
prediction in Otindag Sandy Land 


LUO Jiayan, ZHANG Jing, XU Mengran, MO Yu, TONG Liga 
(1. College of Environment and Resources, Dalian Minzu University, Dalian 116600, Liaoning, China; 


2. College of Architecture, Dalian Minzu University, Dalian 116650, Liaoning, China) 


Abstract: In the context of future global warming, the vegetation in arid and semi-arid regions in China is more 
sensitive, and the balance between economic development and environmental protection remains challenging. 
Therefore, it is essential to conduct dynamic monitoring and driving factors assessment of vegetation in those ar- 
eas. In this study, the spatiotemporal variation of NDVI in Otindag Sandy Land from 2000 to 2020 was first ana- 
lyzed using Theil-Sen median trend analysis; next, the driving mechanism of NDVI change was analyzed using 
the Geodetector model. Finally, multiscenario analysis was used to reveal the future change trajectories of vegeta- 
tion in the area. The results showed that the NDVI of Otindag Sandy Land exhibited upward fluctuation trends 
within 20 years, and grazing is the main factor affecting the changes in NDVI in Otindag Sandy Land; however, 
the effect of grazing to NDVI dynamics in the Otindag Sandy Land gradually weakened. In addition, the multisce- 
nario prediction results show that, although the vegetation of Otindag Sandy Land is developing toward restora- 
tion, its fragile ecological environment and the vegetation status were still threatened by human activities and cli- 
mate change. This study provides a theoretical basis for the Otindag Sandy Land ecological construction and res- 
toration, as well as relevant countermeasures and suggestions according to the actual situation. 


Key words: trend analysis; geodetector; multi-scenario prediction; NDVI; Otindag Sandy Land 


